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ABSTRACT 

We present a coherent timing analysis of the 2003 outburst of the accreting milUsecond pulsar 
XTE J1807— 294 . We find an upper limit for the spin frequency derivative of < 5 x 10~^'*Hz/s. 
The sinusoidal fractional amplitudes of the pulsations are the highest observed among the accreting 
millisecond pulsars and can reach values of up to 27% (2.5-30 keV). The pulse arrival time residuals 
of the fundamental follow a linear anti-correlation with the fractional amplitudes that suggests hot 
spot motion over the surface of the neutron star both in longitude and latitude. An anti-correlation 
between residuals and X-ray flux suggests an influence of accretion rate on pulse phase, and casts 
doubts on the use of standard timing techniques to measure spin frequencies and torques on the 
neutron star. 

Subject headings: stars: individual (XTE J1807— 294 ) — stars: neutron — X-rays: stars 
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1. INTRODUCTION 

An open problem in the field of accreting millisecond 
pulsar (AMXP) is how to devise a reliable method to 
measure spin and orbital parameters. Since the discovery 
of the first AMXP (Wijnands & van der Klis 1998) con- 
siderable improvements have been made, leading to the 
measurement of accurate orbital and spin parameters for 
9 of the 10 known AMXPs (see Wijnands 2004, Poutanen 
2006, di Salvo et al. 2007 for a review and Patruno et al. 
2009 for the last source discovered). Current methods 
(see e.g. Taylor 1992) are based on folding procedures 
to reconstruct the pulse profiles of the accreting neutron 
star and on direct measurement of the pulse phase vari- 
ations due to orbital Doppler shift and spin changes (for 
example due to torques) . The pulse phases are fitted us- 
ing minimization techniques. However, a substantial 
complication sometimes arises due to the presence of a 
strong unmodeled noise component in the pulse phases 
that, when ignored, might affect the reliability of the 
method. Two possible strategies have been used in the 
literature to try and overcome this: (i) harmonic data se- 
lection (Burderi et al. 2006, Riggio et al. 2008, Chou et al. 
2008, Papitto et al. 2007) and (ii) use of a minimum vari- 
ance estimator (Boynton & Deeter 1985, Hartman et al. 
2008). In the first case the pulse profiles are decomposed 
into their harmonic components: generally one sinusoid 
at the fundamental frequency (or first harmonic, v) and 
one at the second harmonic (2z/) and are analyzed sep- 
arately, measuring two independent sets of orbital and 
spin parameters. The harmonic with the weakest noise 
content is selected for the measurement of the spin and 
orbital parameters and the noisier one is discarded (see 
e.g., Burderi et al. 2006). Although this use of the most 
"stable" harmonic reduces the x^j this selection throws 
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away part of the information and in that sense is not op- 
timal. The hypothesis behind the selection of the most 
stable harmonic is that, for unknown reasons, that har- 
monic better tracks the spin of the neutron star. Burderi 
et al. (2006) speculated that the second harmonic might 
be more stable because it arises from accretion onto both 
the polar caps and hence is insensitive to the fiux ratio 
between poles. 

In the second method, both harmonics are used and 
weighted to minimize the effect of phase noise (Boyn- 
ton & Deeter 1985, Hartman et al. 2008). However, 
also in this second situation in practice data selection 
is performed. If the phases of both harmonics change 
differently, the possibility of defining pulse arrival times 
breaks down and the data where this happens have to be 
excluded from the analysis (Hartman et al. 2008). 
Because both methods employ different data selections, 
different results are obtained when analyzing the same 
source. For example in the case of SAX J1808. 4-3658, 
the pulse frequency derivative measured from only the 
second harmonic was 4.4 x 10~"'^^Hzs~^ for the first 14 
days of the 2002 outburst and -7.6 x lO'^'^Hzs-i for 
the rest of the outburst (Burderi et al. 2006). In Hart- 
man et al. (2008) we considered the same source and 
gave an upper limit of \v\ < 2.5 x 10~"'^^Hzs~^ for all 
the four outbursts for which high resolution timing data 
was available. The reason for this discrepancy is that 
while Burderi et al. (2006) used only the information 
carried by the second harmonic and rejected the results 
of the fundamental frequency, we used both harmonics 
but excluded the initial data where the phase variations 
where stronger and discrepant between harmonics (Hart- 
man et al. 2008). So these differences arise as a conse- 
quence of different data selections. In this paper we try 
to better characterize the timing noise such as observed 
in AMXPs focusing on a source where the noise is strong: 
XTE J1807-294 (J1807 from now on) which has been in 
outburst for w 120 days in 2003 (Markwardt et al. 2003). 

2. DATA REDUCTION AND RECONSTRUCTION OF THE 
PULSE PROFILES 

We reduced all the pointed observations from the 
RXTE satellite taken with the Proportional Counter Ar- 



2 



Patruno et al. 



ray (PCA, Jahoda ct al. 2006) that cover the 2003 out- 
burst of J1807. The PCA instrument provides an array 
of five proportional counter units (PCUs) with a collect- 
ing area of 1200 cm^ per unit operating in the 2-60 keV 
range and a field of view with a FWHM of ~ 1°. 

We constructed the X-ray lightcurve using the counts 
in PCA Absolute channels 5-67 (w 2.5 - 29 kcV). 

We constructed our pulse profiles by folding 512 s long 
chunks of lightcurve in profiles of = 32 bins, with the 
ephemeris of Riggio et al. (2008). In this folding process 
we used the TEMPO pulsar timing program to generate 
a scries of polynomial expansions of the ephemeris that 
predict the barycentcred phase of each photon detected. 
The total number of photons detected in a single profile 
bin is Xj ± ^/xj, with the error calculated from counting 
statistics and j = 1, ...,N. Since the pulse profile shape 
changes throughout the outburst, it is not possible to 
base the analysis on a stable template profile. Therefore 
we decided to analyze the pulse profile harmonic compo- 
nents separately. 

To calculate the pulse fractional amplitudes and phases 
we decomposed each profile as: 
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by using standard minimization techniques. The 
term bk is the amplitude of the sinusoid representing the 
fc— th harmonic, and ba is the unpulsed flux component. 
We choose the first peak of each sinusoid in the profile as 
the fiducial point for each harmonic. Defining the fc-th 
harmonic frequency to hek-v, the unique pulse phases (pk 
of each harmonic range from to 1 . The i-th pulse time 
of arrival (TOA) of the k th harmonic is then defined as: 
tk^i = k^ + ^ti- Here Ati is the time of the middle of 
the i-th folded chunk. With these definitions, a positive 
time shift is equivalent to a lagging pulse TOA, while 
a negative shift corresponds to a preceding pulse TOA. 
This is the convention that will be used later to define 
pulse phase residuals. 

The fractional sinusoidal amplitude of the i-th pulse 
profile and the k-th harmonic is calculated as: 
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where ^ph^i and B, are the total number of photons 
and the backgromid counts (calculated with the FTOOL 
pcabackest) in the i-th pulse profile. The error on the 
fractional amplitude Ri^k is calculated propagating the 
errors on bk and Nph^i. The error on Bi is negligible 
with respect to the other errors and will not be consid- 
ered further. 

We define a pulse profile harmonic to be significant 
if the ratio between the amplitude bk and its statistical 
error a^^ is larger than 3.3 when using a folding time 
of 512 s. The choice of 3.3 guarantees that the number 
of false detections expected when considering the global 
number of pulse profiles (» 850), is less than one. The 
length of the folding time was then changed to 300 and 
3000 s to probe different timescales (see §3), and the 
significance threshold rescaled to 3.5 and 3(T respectively, 
according to the new number of pulse profiles. 

After obtaining our set of TOAs for all the signifi- 
cant harmonics we chose to describe the phase (j) of the 



k-th harmonic (we omit the k index from now on) at 
the barycentric reference frame, as a combination of six 
terms: 
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where (t)L (t) is a linear function of the time {(pL (t) = 
+ vt, with (/)o an initial reference phase), (f)Q (t) is a 
parabolic function of time {(pg {t) = ^i^t^), and 0o (0 is 
the keplerian orbital modulation component. The term 
(pM {t) is the measurement error component, and is given 
by a set of independent values and is normally distributed 
with an amplitude that can be predicted by propagating 
the Poisson uncertainties due to counting statistics. The 
term (j)A (t) is the astromctric uncertainty position error, 
and the last term, <j)N (t), is the so-called timing noise 
component that defines all the phase variations that re- 
main. The timing noise includes, but is not limited to, 
any phase residual that can be described as red noise and 
possible extra white noise in addition to that described 
by the measurement error component 0m {t)- 

One of the key points when dealing with timing noise 
is how to distinguish a true spin frequency change of the 
neutron star from an effect that mimics it. In general, 
(j)Q and (f)N- can both be due to torques, both not be 
due to torques, or one can, while the other is not. In the 
first case the torque is not constant and has a fluctuating 
component. In the second case there is a process different 
from a torque affecting the pulse phases. In the third 
case, if (/>q is due to a torque, it is constant, while if (^jv 
is due to a torque then the torque is not constant. 

In the presence of timing noise {(j>N) the formal pa- 
rameter errors estimated using standard minimization 
techniques are not realistic estimates of the true uncer- 
tainties, as the hypothesis behind the minimization 
technique is that the source of noise is white and its am- 
plitude can be predicted from counting statistics. In the 
presence of an additional source of noise, such as the tim- 
ing noise, the apparently significant measurement of a 
parameter can simply reflect the non realistic estimation 
of the parameter errors. To solve this, we adopted the 
technique we already employed in Hartman et al. (2008), 
who used Monte Carlo (MC) simulations of the timing 
residuals to account for the effect of timing noise on the 
parameter errors. The technique uses the power density 
spectrum of the best-fit timing residuals of a ;> model, as 
output by TEMPO. Then thousands of fake power den- 
sity spectra are produced, with Fourier amplitudes iden- 
tical to the original spectrum and with random uncorre- 
lated Fourier phases. The Fourier frequencies are then 
transformed back to the time domain into fake residuals, 
and thousands of v and values are measured to cre- 
ate a Gaussian distribution of spin frequencies and spin 
frequency derivatives. The standard deviations of these 
distributions are the statistical uncertainties on the spin 
frequency and derivative. For a detailed explanation of 
the method we refer to Hartman et al. (2008). 

3. RESULTS 

3.1. Measurement of the spin frequency and its 
derivative in the presence of timing noise 

We fitted the phases of each harmonic with a circu- 
lar keplerian model {(j>o) plus a linear term {(^l) and a 
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quadratic term ((/>q). All the residual phase variation 
we observe after removing these three terms is treated as 
noise (0m and (Pn )- The v and v measured for each of the 
two harmonics is given in Tables 1 and 2, respectively. 
The errors on the pulse frequency and its derivative are 
calculated performing 10* MC simulations as described 
in § 2. At long periods (days), red noise dominates the 
power spectrum, while at short periods (hours), the un- 
correlated Poisson noise dominates. The red noise power 
spectrum is not very steep, and has a power law depen- 
dence P{v) oc I/" with a « —0.5. 

The source position we used comes from Chandra ob- 
servations whose 68% confidence level error circle is 0".4 
in radius. The astrometric uncertainty introduced in 
this way on the frequency and frequency derivative is 
3 X 10^^ Hz and 0.7 x 10"^*Hzs~^, respectively (calcu- 
lated with eqs. Al and A2 from Hartman et al. (2008), 
which added in quadrature to the MC statistical errors 
gives the final errors reported in Tables 1 and 2. The fi- 
nal pulse frequency derivative significances for the funda- 
mental and the second harmonic are ~ 2.7cr and w l.Scr, 
respectively. 

We note that the significance of the frequency deriva- 
tive for the fundamental increases above the 3 sigma 
level when the statistical errors are calculated with stan- 
dard minimization techniques, consistently with Rig- 
gio et al. (2008). These errors calculated with Ax^ = 1.0 
are 2 x lO-^^Hzs-^ and 1.6 x lO^^^Hzs^i for the fun- 
damental and second harmonic respectively. So, a sig- 
nificant v is present which is, however, consistent with 
being part of the (red) timing noise. 

The timing residuals obtained after removing a u — Q 
model are plotted in Figure 1 for both the harmonics (see 
Tables 1 & 2 for the pulse frequencies used in the fits). 
Our orbital solution is consistent for the two harmonics 
and with the orbital parameters published in Riggio et al. 
(2008). For the fundamental we find: 

• orbital period: 2404.4163(3) s 

• projected semi-major axis: 4.830(3) It-nis 

• time of ascending node: MJD 52720.675601(3) 

where the quoted errors are calculated with the min- 
imization technique and correspond to Ax^ = 1. Since 
the pulse phase residuals are approximately white and 
consistent with the expected Poissonian uncertainty, on 
timescales equal to and shorter than the orbital period, 
the orbital parameter errors are a good approximation 
of the true uncertainties. 



3.2. Relation between timing residuals and X-ray flux 

In this section we analyze the relation between the 
pulse arrival time residuals relative to a constant pulse 
frequency {v — 0) model and X-ray flux. Riggio et al. 
(2008) found that the residuals of the fundamental show 
a strong correlation with the X-ray flux, while the second 
harmonic shows only a marginal correlation. Since large 
pulse phase shifts are often observed (in both harmonics) 
in coincidence with the flaring states, we investigate the 
possibility that at least part of the observed timing noise 
is correlated with the presence of X-ray flux variations. 
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Fig. 1. — a): Timing residuals for a constant spin frequency 
and a circular keplerian orbit. The fundamental (blue circles) and 
the second harmonic (red squares) phases were measured using an 
integration time of 512 s per pulse profile. 

b) : Sinusoidal fractional amplitude of the fundamental (blue aster- 
isks: flaring; black circles: non-flaring) and second harmonic (red 
squares) during the whole outburst. During the flares, the funda- 
mental sinusoidal fractional amplitude grows up to 27%, which 
is the highest value ever observed for an AMXP. 

c) : XTE J1807-294 hghtcurve of the 2003 outburst. The count 
rate was normalized to the Crab (Kuulkers et al. 1994) using the 
data nearest in time and in the same PCA gain epoch (e.g., van 
Straaten et al. 2003). The blue circles and the black asterisks 
identify the 4 non-flaring and the 3 flaring states, respectively, as 
deflned in Chou et al. (2008). 

In this section we show that both the harmonics are 
consistent with being correlated with X-ray flux. First 
we focus on the entire data set, then we split the data 
in intervals choosing the same 7 chunks as Chou et al. 
(2008); see Figure Ic), distinguishing non-flaring states 
following the exponential flux decay of the overall out- 
burst, and flaring states, comprising the six spikes in the 
lightcurve. In Figure 2 we plot the residuals vs. the count 
rate for both the fundamental and the second harmonic. 

We applied a Spearman rank correlation test to the flux 
anti-correlation for each harmonic. We accept the null 
hypothesis (no correlation in the data set) if the proba- 
bility p > 1%. If we do not make any data selection, the 
Spearman test shows no correlation in either harmonic. 
However, a clear split in the data is apparent at around 
7 mCrab: below this threshold the residuals seem to fol- 
low a correlation with the flux, while above this thresh- 
old an anti-correlation is visible for both harmonics. The 
Spearman coefficients for the points above the threshold 
are p = —0.8 and p = —0.65 for the fundamental and the 
second harmonic respectively (p < 1%). A few outliers 
are visible in the plot, such as for example the four points 
of the second harmonic at about « —0.3 cycles. These 
points correspond to data taken during some of the flar- 
ing states. If we consider only the non-flaring states, the 
Spearman coefficients become p — —0.9 and p = —0.76 
respectively {p < 1%). 

The fact that we see a change from correlation to anti- 
correlation around 7 mCrabs is due to the fact that at 
that flux level in the decay of the outburst the timing 
residuals reach the peak of the parabolic function that 
dominates the residuals (at MJD« 52745, see Figure 1). 
This is a consequence of the fitting procedure, which se- 
lects the constant reference pulse frequency that mini- 
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TABLE 1 

Timing parameters for XTE J 1807-294 (Fundamental) 



Parameter Fundamental MC error Astrometric error Final error 
Spin frequeney v (Hz) 190.62350702 2 X lO^^'Hz 3 X lO^^'Hz 4 X lO^^Hz 
Spin frequeney derivative i>(10^^*Hzs^^) 2.7 0.7 0.7 1.0 
Reference Epoch (MJD) 52720.0 



TABLE 2 

Timing parameters for XTE J1807-294 (second harmonic) 



Parameter second harmonic MC error Astrometric error Final error 
Spin frequeney v (Hz) 190.62350706 3 X lO^^Hz 3 X lO'^'Hz 4 X lO'^Hz 
Spin frequeney derivative i> (lO^^^Hz s^"^) 1.6 0.8 0.7 1.1 
Reference Epoch (MJD) 52720.0 
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Fig. 2. — Phase residuals vs. X-ray flux for the fundamental 
(blue asterisks for the flaring intervals and open black circles for the 
non-flaring intervals) and the second harmonic (red open squares) 
relative to a i> = model. Each pulse is a 512 s folded chunk of 
lightcurve. The dashed line at around 7 mCrab splits the diagram 
in two regions: in the left one the points follow a correlation, in 
the right one they follow an anti-correlation. The four points of 
the second harmonic that lay outside the relations correspond to 
the large jump observed during the second flare. 

mizes the of the timing residuals. As the observed 
pulse frequency is increasing, the reference frequency is 
too fast for the rising part of the residuals, and too slow 
for the decreasing part. 

We have seen in § 3.1 that the measured pulse fre- 
quency increase is consistent with being part of a red 
noise process and that true neutron star spin variations 
may or may not be the cause. We can choose a higher 
reference pulse frequency than the one used to produce 
Figure la, and turn the correlation-anti-correlation di- 
chotomy in the flux-residual diagram into only an anti- 
correlation, at the cost of increasing by a factor w 10. 
A V higher by 10~^ Hz makes the split in the data disap- 
pear and increases the degree of correlation between flux 
and timing residuals. 

All the correlations and anti-correlations disappear or 
are strongly reduced for the timing residuals relative to 
the best-fit finite constant-i' model. 

3.3. Pulse profiles 



In this section we focus on the shape of the pulse pro- 
files and their relation with other observables, such as 
the phase, the timing noise and the X-ray flux. 

3.3.1. The fractional amplitude-residual diagram 

We have seen in the previous section that for some data 
selections the X-ray flux correlates with the timing resid- 
uals relative to a = model, but not when a finite v is 
admitted. As already noticed by Zhang et al. (2006) and 
Chou et al. (2008), the fractional amplitude of the pulsa- 
tions shows six spikes coincident with the six flares in the 
lightcurve. Therefore, a correlation might also exist be- 
tween the fractional amplitude of the pulsations and the 
arrival time residuals. Using a, ij = model, and again 
using a Spearman rank test, we found a correlation co- 
efficient p = —0.61 {p < 1%) for the fundamental, while 
no significant correlation exists for the second harmonic. 
The second harmonic is also inconsistent with following 
the same correlation as the fundamental. Repeating the 
test for a i/ model we still find no correlations for the 
second harmonic, but the anti-correlation found for the 
fundamental becomes stronger {p — —0.80, p < 1%). In 
Figure 3 we show the fractional amplitude vs. residual 
diagram (relative to a j> model). The anti-correlation is 
evident. It is interesting that the small number of points 
(circled in the figure) that are outliers all belong to the 
first 2.5 days of the outburst. 

We then analyzed the flaring and non-flaring states 
separately. The non-flaring state shows a weak anti- 
correlation with a J> = model {p — —0.43, p < 
1%) which becomes slightly stronger with a model 
{p = —0.51, p < 1%) The flaring state shows an anti- 
correlation relative to a z> = model {p — —0.58, p < 1%) 
that becomes much stronger for a i> model (p = —0.81, 
p < 1%). 

We found no energy dependence in this fractional 
amplitude-timing residual anti-correlation (amplitude 
anti-correlation from now on) when we repeated the anal- 
ysis in 6 different energy bands from 2.5 to 30 keV. The 
same is true for the second harmonic: no correlation was 
found in any energy band. 

3.3.2. The X-ray flux and the fractional amplitude 

In our previous paper (Hartman et al. 2008), we 
found an anti-correlation between the fractional ampli- 
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Fig. 3. — Timing residuals vs. fractional amplitude diagram. 
The blue asterisks refer to the flaring states, while the black circles 
are the non-flaring states, both referring to the fundamental fre- 
quency. The second harmonic is plotted as red open squares. Each 
pulse was built using 512 s of integration time. The residuals are 
relative to a finite v model. The green circled outliers of the anti- 
correlation for the fundamental, all belong to the first 2.5 days of 
the observations. The second harmonic amplitude is uncorrelated 
to timing residuals. 

tude of the second harmonic and the X-ray flux in SAX 
J1808.4— 3658. We also noted that the fractional ampli- 
tude of the fundamental behaved unpredictably. Some- 
thing similar applies to J1807, where no correlation is 
found for the fundamental while a strong anti-correlation 
exists between the observed count rate and the fractional 
amplitude of the second harmonic (p = —0.79, p < 1%, 
see Figure 4). The behavior of the fundamental is in- 
consistent with this relation. By analogy with Hart- 
man et al. (2008) we fitted a simple power-law model 
(i?2 c)c /J, where ]x is the X-ray flux) to the data, 
which gives a power law index 7 = —0.41 ± 0.04 with 
a x^/dof of 90.2/117. Interestingly, the power law index 
we found for SAX J1808.4-3658 (Hartman et al. 2008) 
was in agreement with this. So, a difference in behavior 
exists between the fractional amplitude of the fundamen- 
tal frequency and of the second harmonic. They respond 
differently to both the flux and the arrival time residuals. 



3.3.3. Fractional amplitude 

We focus now on the energy dependence of the pulse 
profiles. We consider again all the data available and the 
subgroups of flaring and non-flaring states. Chou et al. 
(2008) already reported on the energy dependence of the 
fundamental frequency during the non flaring state. Here 
we explore also the flaring state and the energy depen- 
dence of the second harmonic. Looking at Figure 5 two 
interesting features are immediately apparent: 

1. the fractional amplitude energy dependence is the 
same for both harmonics and regardless of the state of 
the source (flaring, non-flaring), up to a constant factor 

2. the fractional amplitude of the fundamental increases 
by a factor of « 1.8 during the flaring state with respect 
to the non-flaring state, while it remains approximately 
constant for the second harmonic. 

Another important property of the pulses is the time 
dependence of the fractional amplitude. In the middle 
panel of Figurel we plot the fractional amplitude of the 



Fig. 4. — The fractional amplitude of the second harmonic is 
anti-correlated with the fiux and scales with a power law of index 
7 = —0.41 it 0.04, close to the power law index found in a similar 
relation for SAX J1808. 4-3658. 

pulsations in the 2.5 — 30 keV band. As can be seen, dur- 
ing the last of the six flaring states the fractional ampli- 
tude of the fundamental increases up to « 27%, which is 
the highest ever observed for an AMXP.'* Selecting a nar- 
rower band between 2.5 and 10 keV the maximum frac- 
tional amplitude does not appreciably change. During 
the non-flaring stage the fractional amplitude decreases 
smoothly from « 9% down to « 4%. The second har- 
monic amplitude on the contrary increases from k, 2% 
up to « 5%. 

3.3.4. Harmonic content 

We decomposed each pulse profile in its harmonic com- 
ponents to look for the presence of higher harmonics. 
While the detection of the second harmonic is quite com- 
mon among the AXMPs, higher harmonics have never 
been detected, with the exception of a possible third har- 
monic in SAX J1808.4-3658 (Hartman et al. 2008). In 
J 1807 we detected a third harmonic at better than 3a, 
in several different stages of the outburst, with a max- 
imum fractional amplitude of w 1.5% at MJD around 
52560. To increase the S/N, we folded chunks of data 
of length 3000 s. The number of > 3(7 detections of the 
third harmonic was of 11 out of 163 chunks searched. 
We searched the same chunks for a fourth harmonic, and 
found 5-10 significant detections above 3ct in the whole 
outburst, depending on the binning. When detected, the 
fourth harmonic has a fractional amplitude 0.5 — 2.0%. 

There were no observations where we detected all 4 
harmonics at the same time. During the second and third 
flares, we found a second and fourth but not a third 
harmonic, during the first two flares we found a second 
and third but not a fourth harmonic. 

For the third and fourth harmonics we count respec- 
tively 8 and 5 detections during the flaring states and 3 
and 2 detections in the non-flaring states. 

The fractional amplitude of the third harmonic also de- 
creases with the flux, although the slope of the power law 
is much smaller (7 = -0.017 ± 0.004). The fourth har- 
monic has no significant fiux dependence, but its power 

^ In this paper we are quoting sinusoidal fractional amplitudes, 
which are \/2 larger than the rms fractional amplitudes 
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law slope is also consistent with the 7 obtained for the 
third harmonic. 

Of course this result has to be taken with caution, 
since we are suffering from low number statistics with 
only « 20 detections of the third and fourth harmonic 
altogether. 
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Fig. 5. Energy dependence of the pulse fraetional amplitudes. 
The squares and the triangles refer to the flaring and non-flaring 
states respectively. The circles comprise the whole outburst. The 
bottom curves, overlapped in the plot, are the fractional amplitudes 
of the second harmonic which remains stable in both states. The 
pulses of the fundamental in the flaring states have a fractional 
amplitude which is about 1.8 times larger than during the non- 
flaring states. Up to a constant factor, the fractional amplitude 
has the same energy dependence for both harmonics and for both 
flaring and non-flaring states. 



3.4. Short-term v measurements 

Using the fundamental frequency, we measured short- 
term pulse frequency derivatives using the seven sub 
groups of data as defined in §3.2. These measurements 
are useful to investigate the time dependence of the pulse 
frequency derivative with time. This test is possible in 
J1807 because of its very long outburst duration (more 
than 120 days, of which « 106 days with detectable pul- 
sations). 

The V values and their uncertainties were first calcu- 
lated with standard minimization techniques. All 
measured v values during the non-flaring states had a 
positive sign, whereas a negative sign was measured for 
all three flaring states. The measured v values are shown 
in Figure 6. There is no clear trend, and most impor- 
tantly no correlation between ii and the average X-ray 
flux in either the flaring and the non-flaring states. This 
test cannot be performed on the second harmonic, since 
the smaller number of detections prevents a meaningful 
analysis of data subsets for this purpose. 

We then calculated the statistical uncertainties on the 
i> for each sub group of data by using the MC method as 
explained in § 3.1. All the i/ values were consistent with 
being part of the same red noise process, consistently 
with what was calculated for the long term v value of 
§ 3.1. 

4. DISCUSSION 

We have analyzed the outburst of XTE J1807-294 and 
we have calculated statistical errors by means of MC 
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Fig. 6. — Frequency derivative (i>) evolution. The non flaring 
states (open circles) all have positive !>, that however does not fol- 
low a power-law decrease as expected from the accretion theory. 
The flaring states (asterisks) all have a negative value, correspond- 
ing to spin down. 

simulations as we previously did for SAX J1808. 4-3658 
(Hartman ct al. 2008). We found that with our statis- 
tical treatment of the red noise observed in the timing 
residuals of both the fundamental and the second har- 
monic, the significance of the spin up is reduced below 
3(7 for both the fundamental and the second harmonic. 

The fact that the spin frequency derivative is not sig- 
nificant does not mean that there is not a component in 
the residuals that can be fitted with a parabola. It just 
means that the parabola is consistent with having the 
same origin as the power at other low frequencies: both 
the parabola and the remaining fluctuations are consis- 
tent with being part of the realization of the same red 
noise process in the timing residuals. It is a separate issue 
whether or not this process is due to true spin changes 
and torques on the neutron star. 

Our observed parabola in the timing residuals com- 
bined with the stochastic and astrometric uncertainty 
implies that any spin frequency derivative has a magni- 
tude smaller than |z/| ;$5 x 10~^^Hzs~^ at the 95% con- 
fidence level. 

Evidence against the spin-up interpretation of the 
phases comes from the lack of any correlation between 
the observed X-ray flux and the measured v (see § 3.4). If 
standard accretion torque theory applies, then the mag- 
netospheric radius (r^) should decrease as the mass ac- 
cretion rate M increases, following a power-law oc 
M ^" when r,„ < r-co, with a = 2/7 in the simplest case, 
where r^o ls the corotation radius. This implies that also 
the instantaneous v has a power-law dependence on the 
mass accretion rate, and (when < fco) it is: 



(4) 



see Bildsten et al. (1997). Here M is the average mass 

accretion rate, M the neutron star mass and / the neu- 
tron star moment of inertia. We have observed no such a 
correlation between the flux and the instantaneous pulse 
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frequency derivative, neither in the flaring nor in the non- 
flaring states. One possible explanation is that the X-ray 
flux is not a good tracer of the mass accretion rate. If 
it is, standard accretion theory does not apply and the 
most logical conclusion is that the observed timing resid- 
uals are not due to torques. 

The possibility that the X-ray flux is not a good tracer 
of the mass accretion rate is a long standing issue in 
the X-ray binary pulsar field and has no simple solution. 
If the X-ray flux is completely unrelated to the mass 
accretion rate, then no conclusions can be drawn on the 
eff'ect of the accretion on the pulse phase. 

By using eq. (4) we can calculate the spin frequency 
derivative expected for J1807 from standard accretion 
theory, assuming a distance of 8 kpc and converting the 
average X-ray luminosity into an average mass transfer 
rate through Lx ~ r)(?M. We assume; an efficiency q = 
0.1 for the conversion of gravitational potential energy 
into radiation. In this way we obtain an average mass 
accretion rate M « 3 x 1O~"'^'^M0 yr""'^ (averaging over 
the outburst). Assuming « rco we have an expected 
z> « 10~^^Hzs^^, which is below our calculated upper 
limit of 5 X 10~^''Hzs~^. However, the short term 
values calculated in §3.4 exceed the theory value by 1-2 
orders of magnitude and therefore are very unlikely due 
to accretion torques. 

The possibility that we are not observing the effect of 
a torque on the neutron star is also suggested by the fact 
that looking at the shape of the lightcurve one can im- 
mediately infer the sign of the measured pulse frequency 
derivative in the timing residuals. This is a consequence 
of the fiux anti-correlation. If the lightcurve is concave, 
then the average i> is positive, while if the lightcurve has 
a convex shape the average will be negative. This ex- 
plains why ;> > in the non-fiaring states and i) < 
in the flaring states. It suggests a direct influence of 
the accretion rate on phase, which could be effectuated 
through the hot spot position on the neutron star sur- 
face. Extending this interpretation to the average i/ over 
the entire outburst, we also favor the interpretation of a 
moving hot spot for that long term trend, discarding the 
hypothesis of a torque to explain the parabola observed 
in the pulse phase residuals. 

Chou et al. (2008) also suggested that the lagging ar- 
rival times observed during the flaring states cannot be 
explained with a torque model, since they correspond to 
a sudden change from a spin up to a spin down. These 
authors also suggested that motion of the hot spot can be 
responsible for both the phase shifts and the increase of 
the fractional amplitude during the flaring states. Chou 
et al. (2008) assumed a fixed position of the hot spot 
during the non-flaring states. However, it is unlikely that 
the hot spot is fixed on the surface during the non-fiaring 
state, as we have shown (see § 3.4) that the magnitude 
of the short-term i/ is too large to be compatible with 
standard accretion theory. 

Ibragimov & Poutancn (2009) recently proposed a re- 
ceding disk as a possible explanation for the timing noise 
and pulse profile variability observed in the 2002 outburst 
of SAX J1808.4-3658. In this model the antipodal spot 
can be observed when the inner accretion disk moves 
sufficiently far from the neutron star surface as a conse- 
quence of decreasing flux. We observed a strong overtone 



and pulse phase drifts since the early stages of the out- 
burst, when the disk should be closest to the neutron 
star. Therefore it is not clear whether our observations 
can be explained by this model or not, and further in- 
vestigations of the problem are required. 

Two hot spots with different and variable intensities 
can produce a phase shift and a changing pulsed ampli- 
tude, even if the location of both hot spots is fixed on 
the neutron star surface (Burderi 2008). This possibil- 
ity needs also further investigation since a self- consistent 
model has not yet been presented. 

We observed (1) a relation between flux and time of 
arrivals for both the flaring and non-flaring state (§ 3.2). 
This relation was consistent with being the same for 
the two states. We also observed (2) an anti-correlation 
between pulse fractional amplitudes and time of arrivals 
during the flaring state. Finally, (3) this amplitude 
anti-correlation became stronger when using a long 
term z>. The amplitude anti-correlation was weak in the 
non-flaring state, regardless of the i>. In the context of a 
hot spot motion model for the time of arrival variations, 
these findings constrain the kinematics of this motion. 

Lamb et al. (2008) demonstrated that variations in 
the pulse fractional amplitudes should be anti-correlated 
with their time of arrival if the hot spot is close to the 
neutron star spin axis and the hot spot wanders by a 
small amount in latitude. 

Lamb et al. (2008) showed that even a small displace- 
ment in longitude of the emitting region, when close to 
the spin axis, produces a large phase change, but no am- 
plitude variation. A motion in latitude produces both 
phase and amplitude changes due to the hot spot veloc- 
ity variation affecting Doppler boosting and aberration. 
An anti-correlation between the pulse arrival times and 
the pulse amplitudes would be an indicator of the above. 
Combining this with our observational findings 1-3 above 
we conclude within the moving hot spot model for the 
phase variation that: 

1. the hot spot moves with flux in both flaring and non- 
flaring state, since the relation between flux and arrival 
times is observed in both cases and is consistent with 
being the same, 

2. the amplitude anti-correlation in the flaring state 
implies an hot spot moving in latitude. The hot spot 
cannot move mainly in latitute during the non-flaring 
state since a weak amplitude anti-correlation is observed 
and the fractional amplitude changes by only a factor 
« 2 in 106 days. 

3. The long term u must be related with a motion 
in longitude since during the flaring state the amplitude 
anti-correlation becomes much stronger when a i> model 
is used to fit the time of arrivals. This is also compati- 
ble with the non-fiaring state, since the amplitude anti- 
correlation remains weak with or without a 0. 

4. Finally a motion in longitude during the flaring state 
or in latitude during the non-flaring state is possible, but 
it has to be small enough to preserve the observed flux 
and amplitude anti-correlations. 

The reason why the hot spot should drift mainly in 
longitude during the non-flaring states and mainly in 
latidude during the flares might be related with differ- 
ences in the accretion flow process. A hot spot motion 
has been observed in MHD simulations, with a compli- 
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cated dependence of the hot spot position on the mis- 
ahgnment angle between magnetic field and rotation axis 
(Romanova et al. 2002, 2003, 2004). As noted by Lamb 
et al. (2008), long term wandering of the hot spot can 
be related with the structure in the inner part of the 
accretion disk and therefore should track the long-term 
changes of the accretion rate. The position of the hot 
spot on the neutron star surface is expected to change 
rapidly and irregularly as the accretion flow from the in- 
ner region of the accretion disk varies. Further studies 
are required to couple our inferred hot spot kinematics 
to physics and geometry of the accretion flow. We note 
that the fractional amplitudes can also change according 
to the hot spot angular size and/or to the difference in 
temperature between the hot spot and the neutron star 
surface. 

The maximum observed sinusoidal fractional ampli- 
tude (Fig. 1: « 27%) can be explained if the hot spot is 
slightly misaligned from the spin axis (colatitude ^20°) 
with an inclination of the observer larger than w 45°, or 
if the inclination of the observer is smaller than « 45° 
but the spot has a large colatitude (see Figure 2 in Lamb 
et al. 2008, note that we quote sinusoidal amplitudes 
while they use rms amplitudes). 

J1807 shows an anti-correlation between the second 
harmonic fractional amplitude and the X-ray flux. We 
observed a similar anti-correlation in SAX J1808.4— 3658 
(Hartman et al. 2008). This suggests the same process as 
the origin of the anti-correlation in both pulsars. Hart- 
man et al. (2008) found that the anti-correlation was a 
signature of the increasing asymmetry of the pulse pro- 
files toward the end of the outburst. In J1807 the second 
harmonic is less often detected in these late stages of the 
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outburst. However, the lower cotmt rates late in the out- 
burst lead to upper limits on the second harmonic there 
that are sufficiently high that the explanation we pro- 
posed for J1808 (Hartman et al. 2008) can still be valid 
for J1807 as well. 

5. CONCLUSIONS 

In this paper we analyzed the 2003 outburst of XTE 
J1807-294 and found that the pulse frequency deriva- 
tive previously reported in literature is consistent with 

being part of a red noise process. No significant spin fre- 
quency derivative is detected when considering this red 
timing noise as a source of uncertainty in the calcula- 
tionn of statistical uncertainties, and an upper limit of 
5 X 10~^^Hzs~^ can then be set for any spin frequency 
derivative. The average accretion torque expected from 
standard accretion theory predicts a long-term spin fre- 
quency derivative which is still compatible with the de- 
rived upper limit and cannot therefore be excluded from 
current observations. 

We propose hot spot motion on the neutron star sur- 
face as a simpler model able to explain all the observa- 
tions reported in this work, as well as the presence of 
a pulse frequency derivative. If this explanation is cor- 
rect, similar flux and amplitude anti-correlations should 
be observed in other AMXPs. 
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